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Abstract 

In our recent paper we described relationships between integrable systems in- 
spired by the AGT conjecture. On the gauge theory side an integrable spin chain 
naturally emerges while on the conformal field theory side one obtains some special 
reduced Gaudin model. Two types of integrable systems were shown to be related 
by the spectral duality. In this paper we extend the spectral duality to the case 
of higher spin chains. It is proved that the iV-site GL/j Heisenberg chain is dual 
to the special reduced k + 2-points gljy Gaudin model. Moreover, we construct an 
explicit Poisson map between the models at the classical level by performing the 
Dirac reduction procedure and applying the AHH duality transformation. 
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1 Introduction 



This paper is a continuation of pQ where the AGT correspondence [2] was studied at the level 
of integrable systems [3j |3J [5] (see also [6]- [22]). Two sides of the AGT relation correspond to 
a priori different types of integrable models which should actually coincide due to the AGT 
correspondence. This leads to non-trivial predictions of equivalence of different models and also 
illuminates what the equivalence exactly means. The full AGT correspondence associates the 
conformal block of the Virasoro or TV-algebra in two-dimensional conformal field theory with 
the LMNS integral [23] (Nekrasov functions [23])) describing the two-parametric deformation 
of Seiberg-Witten theory by f2-background. Classical integrable systems emerge when the both 
deformation parameters are brought to zero, while when only one of the parameters goes to 
zero (the Nekrasov-Shatashvili limit [4]) the integrable system gets quantized [5j. We study the 
correspondence between the AGT inspired integrable systems only in these two limiting cases. 

In pQ the simplest example of this kind was considered: the equivalence of the four-point 
conformal block and the prepotential in the SU (TV) SUSY theory with 2N fundamental matter 
hypermultiplets and vanishing /^-function. On the gauge theory side the (classical) integrable 
system is known [25J to be the Heisenberg chain [26J which is described by the spectral curve 
r Hciscn (w, x) : det(iu — T{x)) = with GL,2-valued iV-site transfer-matrix T{x) and the Seiberg- 
Witten [27] (SW) differential dS Hciscn (w, x) = x^. On the CFT side the corresponding inte- 
grable system was argued to be some special reduced Gaudin model [28] defined by its spectral 
curve r Gaudm (y, z ) : det(y — L(z)) = with gl^-valued Lax matrix L(z) and the SW differential 
dS Gaudin (y,z) = ydz. 

The spectral duality [29\ 130] generalizes the well-known self-duality of the Toda chain [31J, 
[25j and establishes relation between the two types of models in terms of the bispectral involution 
|32j which interchanges the eigenvalue variable and the spectral parameter. In our case the 
spectral duality at the classical level states that the change of variables z = w, y = x/w relates 
the curves and the SW differentials of the two integrable systems: 

pGaudin^ ^ = pHeisen^ ^ ^ 

(1.1) 

dS" Gaudin (y, z) = dS Hciscn (w, x) . 

The quantum version of the duality emerges from the exact quasi-classical quantization of the 
spectral curves based on the corresponding SW differentials: 

f HciBd^ hzd z )^ Hciscn {z) = , (1.2) 

f Gaudin (^fid z ,z)^ GaUdin {z) = (1.3) 

with some choice of ordering. In [1] the quantum spectral duality was proved in the form: 

fHeisen^ feQJ ^ f Gaudin^^ ^ ( L4 ) 

thus, 

ypHeisen^ _ ^Gaudin _ (1-5) 

The purpose of the paper is twofold: first, to extend the results of [1] to the case of 
higher spin chains. It is proved that the iV-site GL^ Heisenberg chain is dual to the special 
reduced k + 2-points gLy Gaudin model (Theorem 2, Section [5]); and second, to construct 
an explicit Poisson map between the models, i.e. to find an explicit change of variables at 
the classical level. For the second purpose, we perform the Poisson reduction procedure via 
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the Dirac brackets (33^ [M] in the Gaudin model and show that the reduced model exhibits 
the quadratic Poisson structure (Propositions 1,2 in Section [3]). Then, we apply the AHH 
duality transformation [29] and prove that the quadratic Poisson algebra of the dual reduced 
Gaudin model coincides with a natural quadratic algebra underlying the Heisenberg spin chain 
(Theorem 3, Section [5|). These results establish the exact equivalence of the models: 



AHH|zL Gaudin (z)}(x) = T l 
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(1.6) 



Here L Gaudin (z) is the Lax matrix of the Gaudin model and T Hoiscn (x) is the transfer-matrix of 
the Heisenberg chain. 

The paper is organized as follows: in Section 2 we review basic definitions and constructions 
related to the spectral dualities in integrable systems. In Sections 3 and 4 the Gaudin model 
and the Heisenberg chain are described in detail. In Section 5 the classical spectral duality is 
proved and the explicit Poisson map between the models is presented. The quantum version of 
the duality is described in Section 6. In the end we discuss some open problems and comment 
on relation of our results to the previously known. 
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2 Spectral Dualities and Integrable Systems 
2.1 Spectral curves and Poisson structures 

Consider a classical algebraically integrable system [35] with g degrees of freedom. We assume 
that it is described by the Lax matrix L(z) G Mat(N) with the spectral parameter z, local 
coordinate on some complex curve X. It means that the equations of motion with respect 
to any of the Poisson commuting Hamiltonians {H a ,Hp} = 0, a, (3 = l...g (Liouville-Arnold 
integrability |36j) have the Lax form: 

d ta L{z) = {H a ,L(z)} = [L(z),M a (z)}. (2.1) 

The Lax equations may be interpreted as the compatibility condition of the linear problem: 

f (L(z) - X)</> = 0, 

\ (d ta +M a (z))<p = o (Z - Z) 

The first equation gives rise to the spectral curve (a ramified A r -fold covering of X) 

T(A,z) : det(A-L(z)) =0. (2.3) 

which encodes the Hamiltonians and coupling constant:^] in the sense that tr L k (z), 1 < k < N 
are generating functions of them. The spectral curve itself does not fix the integrable system. 

1 The coupling constants are the Casimir functions of the corresponding Poisson brackets generated 
by some underlying classical r-matrix structure. 
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Indeed, V itself does not contain any information about the Poisson structure. Moreover, any 
Lax matrix L(z) is denned up to multiplication by an arbitrary function of z and of integrals 
of motion. Therefore, some more ingredients should be added to specify the integrable model. 
A straightforward way to do this is to introduce some classical r-matrix |31 1. [371 [38| [39] which 
defines the Poisson brackets between any matrix elements of the Lax matrix. The basic examples 
of the r-matrix structures are given by the linear 

{L(z) <f L(z')} = \L(z) <g> 1 + 1 <g> L(z), r(z, z')] (2.4) 

and quadratic 

{L(z) f L(z')} = [L(z) ® L(z'), r(z, z')) (2.5) 

brackets. Each of these structures guarantees the involution {Hi^, Hj ^} = for the z-expansion 
coefficients Hik of tr L k (z). 

A general construction (see, e.g., |40| ) provides solutions for the integrable model in terms of 
the theta-functions on the Jacobian variety Jac(r): the Liouville torus of the integrable system^. 
This construction implies, of course, some Poisson (or symplectic) structure and some choice 
of Darboux coordinates on the phase space. In fact, in |40[ [3TT [4"T] the canonical variables are 
chosen in accordance with the separation of variables (SoV) recipe [42j[43]. This recipe provides 
a set of variables with the Poisson brackets 

{X a ,zp} = h a (\ a ,z a )5 a i3, {X a , Xp} = {z a , zp} = 0, a, f3 = l...g, (2.6) 

where h a are some functions. Each pair (X a ,z a ) lies on the spectral curve, i.e. 

T(X a ,z a ) = 0, a = l...g. (2.7) 

Remark: (|2.7p can be viewed as "separated equations" which appear in the SoV method when 

9 

the generating function of the canonical map is taken in the form S = ^2 S a (z a , c±, c g ), 

a=l 

where c Q are fixed values of the Hamiltonians H a that correspond to the given point of the 
moduli space of (spectral) curves. Then equations (12. 7p are equivalent to \ a = By the 

Liouville-Arnold definition, the integrable model is the Lagrangian bundle which base is the 
set of values of the first integrals H a and the fiber is the Lagrangian submanifold isomorphic 
to Jac(r). The existence of the separated variables means that the Lagrangian submanifold 
is globally presented as a product of g copies T x ... x T of the spectral curve. The recipe 
[4"2l |4"3] provides a concrete set of separated variables. These are the poles z Q of the "properly 
normalized" Baker- Akhiezer function <j) (]2.2h and the dual variables (in the simplest case of the 
normalization) are the corresponding eigenvalues A a for (|2,4[) or log \ a for (|2.5p . For example, 
in sl2 case the variables are defined as zeros of £12(2) while the dual ones are the values of 
L n (z a ). 

Therefore, the spectral curve is written in terms of separated variables. Then there rea- 
sonably exists on r a generating differential dS (of "pdq" type) which provides the action 
variables: 

a a = § dS , 

9 Tsw= fdS, a = l,...g. (2 - 8) 



da a 

B, 



where A4 , Bi are suitable cycles (homology basis) over V and J-sw is the prepotential. This dif- 
ferential (which is simply Adz for some cases) appears naturally in the context of Seiberg-Witten 
theory [27] and is often called the Seiberg-Witten (SW) differential. Choosing a particular dS 



2 In concrete cases some additional factorization may be required to make genus(r) = g. 
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one fixes the functions h a in ()2.6p . For example, if dS = \dz then h a = 1. The Poisson structure 
()2.6p can be also given in terms of the holomorphic symplectic form [41] : 



w = <y(dS), (2.9) 

where 5 denotes the exterior differential on the total space of the fibration of spectral curves 
and divisors [z\, z g \. 

The Lax equations (12. II) are invariant with respect to the gauge transformations 



L(z) -> g(z)L{z)g- 1 {z), M a {z) -> g{z)M a {z)g-\z) + d ta g(z)g-\z) . (2.10) 

The corresponding r-matrix of course changes, while the spectral curve and the SW differential 
remain intact. In fact, we deal with a special (still, a wide) class of Lax matrices having 
only simple poles on S. The gauge transformations may change residues and/or produce new 
poles. Then they connect different phase spaces equipped with different Poisson (and r-matrix) 
structures 



Taking all the aforesaid into consideration, one may conclude that there are two ways to 
define the Poisson structure of an integrable model. The first (gauge invariant) one is to choose 
the SW differential, the second (gauge dependent) way is to define an r-matrix structure. We 
will use the first way and together with the spectral curve this defines the integrable system in 
terms of the separated variables. 

Let us briefly illustrate the construction for the SW theory which will be our primary 
focus in this paper. It is the SU(iV) gauge theory with Nf = 2N matter hypermultiplets 
in the fundamental representation. The corresponding classical integrable system is the GL2 
Heisenberg XXX spin chain with N cites [451 [46] and the spectral curve is written as follows: 

!— : P W - J^ K+{x) - ^if-fx) = 0, (2.11) 



where 

P(x) = H(x-d> k ), K± = J] (z - rrit ) ■ (2.12) 



N N 



The SW differential 



k=l k=l 

dS = ^. (2.13) 

w 

In the gauge theory the parameters cftk are the (diagonal) vacuum expectation values of the 
adjoint scalar field appearing in the M = 2 vector multiplet (Y^k^k = 0), m fc denote the 
masses of the hypermultiplets. The effective low-energy prepotential -Fsw is defined by the A- 
and i?-periods of the SW differential f)2. 13j) on the Riemann surface (|2. 1 1 [) : 

/x dw 
, 2.14 
w 

Ai 

dJ^sw f xdw 



dai 

B. 



(2.15) 

J w 



Remark: Notice that the SW differential can be also chosen as dS = — log(w)dx. This differ- 
ential gives the same answer as (|2,13p for the integrals (|2.8p since both differentials correspond 
to the same representative of the cohomology class if 1 (r,C). 
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2.2 Quantization 



There are two natural ways to quantize the integrable model. In accordance with its Poisson 
(symplectic) structure, one can quantize either the r-matrix structure or the Poisson structure 
in separated variables (|2.6p corresponding to the given SW differential. Let us start from the 
second possibility. 

Considering the SW differential as a symplectic 1-form [41] on C 2 -plane (y, z) yields a pair 
of canonical variables (p(y , z) , q{z)) which brings the SW differential to dS(y, z) = pdq. Then 
there is a natural quantization of the spectral curve defined by the rule (p, q) — > (hd q ,q). 
Therefore, the quantization follows from the recipe based on the given SW differential in the 
quasiclassical form: 

with dS = Adz: (A, z) (-ihd z ,z). (2.16) 

with some choice of ordering. This choice may provide /i-corrections to the coefficients of the 
^-expansion of the quantum spectral curve. In the case of (|2.13p . one has: 

T(w,x) = ^ f = %WJ = 0. (2.17) 



The wave functions can be written in terms of the quantum deformation of the SW differ- 
ential on the spectral curve, that is, ^(z) = exp (—1 f q dS(h)), where dS(h) = p(q,h)dq and 
p(q, 0) = p(q)\r- The monodromies of the wave function around A- and B- cycles of T are given 
by the quantum deformed action type variables [5]: 

*(* + Ai) = exp (-±af) ¥(*), a\ = § dS{h) , 

At 

(2.18) 

^(z + J B i ) = exp(-i&)^(z), ^ = §dS{h), 

where J- NS is the Nekrasov-Shatashvili limit [I] of the LMNS integral [23]. It should be men- 
tioned that we do not impose any boundary conditions which provide a valuable quantum 
problem, i.e. we do not specify the wave functions explicitly. Instead, we analyze the differen- 
tial operator of the quantum spectral curve. 

The differential equation in the r.h.s. of (I2.17P is the Baxter equation [17] . One also 
may choose another quantization. For example, (w,x) — > (e~ h9x ,x), i.e. w maps to the shift 
operator. Then, the Baxter equation is written in the difference form (Fourier dual). From 
the point of view of (|2.16p . the latter case corresponds to dS = — log(w)dx. However, this 
differential gives the same answer as (]2. 13|) for integrals (|2.8|) since both differentials correspond 
to the same representative of the cohomology class i? 1 (r,C) as it was mentioned before. 

Originally, the Baxter equation arises within the Quantum Inverse Scattering Method 
(QISM) [HI [39]. The QISM provides quantization of the phase space and the correspond- 
ing Poisson structures (|2.4j) - (|2.5p via 



L\{z)Ll(w)R\ 2 {z,w) = RUz^)L h 2 {w)L\{z), 
R* 2 (z, w) = 1 <g> 1 + hr{z, w) + 0{h 2 ) 



relations, where L\ = L <g> 1, L2 = 1 <8> L. Further development of the QISM requires some 
quantization conditions, the Bethe equations. In our approach we deal with exact quasi-classical 
equations and do not discuss Bethe-like equations because we do not analyze concrete solutions. 

Besides the approach proposed here, different quantizations of the Gaudin model are known 
(see, for example, [49] and [50]). In [50] there was suggested a method for evaluation of det(c^ + 
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L{z)) as the generating function of commuting quantum Hamiltonians. However, this method 
is based on the linear commutation relation in the corresponding Lie algebra. In the classical 
case this corresponds to the linear Poisson-Lie structure. In our case (see below) we deal with 
the reduced Gaudin model which is described by the quadratic brackets. Therefore, the method 
of [50] is non- applicable in our case (or, at least, requires some verification). In this paper we 
use the recipe (|2.16p which provides the Baxter equation, i.e. a natural quantization of the 
spectral curve based on the separated variables. 

2.3 Bispectral problem and p-q duality 

The notion of bispectral differential operators appeared in the works of F.A.Griinbaum [51] and 
J.J.Duistermaat [52J. G.Wilson formulated it as a bispectral problem [32J: construct the linear 
ordinary differential operator L = X]j=o Lj(z)di with a nonempty family of eigenfunctions 
fy(z, A) depending smoothly on the spectral parameter z such that they are also eigenfunctions 
of a linear ordinary differential operator T = YlT=oTr{X)d^ w hh an eigenvalue g(z) which is a 
function of z: 

L(z,d g )*(z,\) = f(\Mz,\), , . 

C(X,d x )^(z,X)=g(z)^(z,X) { • ) 

It appeared that for the Schrodinger operator L = d 2 z + V(z) the simplest solutions to the 
problem are given by V(z) = (Bessel) and V(z) = z (Airy) cases [52]. Less trivial solutions 
can be obtained by applying the rational Darboux transformations. They satisfy the KdV 
equation [53] . After a link to the Calogero type systems was also found [MIES] it became clear 
that the bispectral problem was closely related to the theory of integrable systems [56]. The 
bispectral problem resembles the quantum version of the p-q duality [571 EH] while the case of 
our interest is somewhat different. Indeed, the p-q duality changes the coordinates of the model 
to the action variables of the dual model while the spectral duality exchanges coordinates and 
momenta in separated variables. One can expect a certain relation between these two types of 
dualities since the separation of variables (|2.6h is "close" to the construction of the action-angle 
variables. Indeed, after the variables are separated, the map to the action-angle variables is 
quite simple because it can be made separately for each degree of freedom. 

At the same time, the p-q duality is very different from the spectral one. While the archety- 
pal example of the spectral self-duality is the Toda chain (see this example in Section 2.4), 
the p-q self-dual model is the rational Calogero-Moser system (and also the trigonometric 
Ruisenaars-Sneider and the hypothetic Double Elliptic Model [57]). Moreover, from the group 
theory interpretation of the p-q duality it follows that the dual models possess Lax representa- 
tions of the same size while in the spectral duality they are different (say, 2x2 and iV x N). 

2.4 Spectral duality 

The duality we investigate in this paper is generated by the bispectral involution [32] which is 
simply a change of arguments of the function 

6 Wilson : T(A,z) -> T(z,A) (2.21) 

corresponding to some spectral problem. 

Definition 1 Let a pair of (algebraically) integrable models be described by the spectral curves 
r(A, z) = 0, T'(X',z') = and the corresponding SW differentials dS(X,z), dS'(X',z'). Then 
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the models are called spectrally dual at the classical level if there exists a change of variables 

X' = X'(X,z), z' = z'(X,z) 

such that 



and 



T(A, z) = b mison [V] (A'(A, z), z'(X, z)) = T'(z'(X, z), A'(A, z)) 



dS(X,z) ^ dS{X'{X,z),z'{X,z)), 



(2.22) 



(2.23) 



where = emphasizes that the SW differential for the integrable system is determined up to a full 
differential on the spectral curve. 

Let us give the very well-know 

Example |31J: The periodic Toda chain can be described by both the gl(iV)-valued Lax matrix 



/ Pi 

e 5(92-91) 



jToda ( \ 
^NxN\ z ) 



e |(93 -92) 



|(92-9l) 

e |(93-92) 

P3 



P2 



\ I e 5(91"9Jv) 

and the GL(2)-valued transfer-matrix [31] 

T 2 T x °f (A) = L N (X)...Lt(X), U(X) = 







X — pi e qi 
-e _ « 



ze \(qi-iN) \ 




PN J 



(2.24) 



i = h 



,N 



(2.25) 



The spectral curves defined by these representations are related by the bispectral involution, 
i.e. 

det(A - L{z)) = and det(z - T(A)) = (2.26) 

coincide. The SW differential is the same in both cases d£ = A^p. Therefore, the periodic Toda 
chain is a self-dual model. 

In quantum case we use the quantization scheme (|2.8p with some choice of ordering. 

Definition 2 Let two integrable models be described by the Baxter equations 

f* = and f'*' = 0. (2.27) 
They are called spectrally dual at the quantum level if their Baxter equations coincide. 

In this paper we prove that the special reduced gl N Gaudin model is spectrally dual to the 
XXX Heisenberg chain at the classical and quantum levels. Moreover, we present an explicit 
Poisson map between the models at the classical level. 

Remark: At the classical level, the coincidence of the spectral curves was mentioned in 
|15j for N = 2. For arbitrary N the general form of the spectral curve for the Gaudin model 
was given in [13 . In quantum case the Baxter equation for the gl 2 Gaudin model was derived 
in [12]. 
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3 Gaudin Model 



Let z be a local coordinate on CP 1 . The Lax matrix is a gl^-valued function L G (z) on 
CP 1 \{^i, . . . , z n } with only simple poles at {z±, . . . , z n } and given residues Kes Zc L(z) = A c G 

n AC 
c=l 

The spectral curve is 

f Gaudin (y, z) : det(y - L G (z)) = 0. (3.2) 

3.1 Unreduced Gaudin model 

The phase space of the Gaudin model [28], [59] is a direct product of orbits of the coadjoint 
action of GL^r: 

^Gaudin = d X ... X O n , 



it 



dimWl Gaudin = £ dimO c 

c=l 



(3.3) 



This phase space is equipped with the Poisson-Lie brackets: 

{A^A c kl } = 5 bc {A c kj 8 ll -A^8 kj ) , 6,c=l,...,n, (3.4) 
which is generated by the linear r-matrix structure: 

{L G (z) f L G {w)} = [L G {z) §1 + 1® L G (w),r(z, w)], 



N 



(3.5) 



The orbits are realized by fixation of the Casimir functions or eigenvalues of {^4 C }, i.e. 

A c = gAlg- 1 , g G GLjv , A c = diag(Ai, X N ) e gl* N . (3.6) 

The spectrum j4q defines the dimension of O c . For example, in general case (when Aj are 
arbitrary) 

dimO max = iV(iV-l). (3.7) 

In the case when N — 1 eigenvalues coincide 

dimO min = 2(N - 1). (3.8) 

The later orbit can be parameterized in a "quiver-like" way [60| using the vector (column) £ 
and the covector (row) T] T : 

A = ixr] T , Aij = fay, {&, rij} = Sij (3.9) 

The symmetry 

^ a£ v .+ i q (3.10) 

a 

N 

generates the "conservation law" ^ 77^ = tiA = NX = const. Fixing the gauge as tjn = 1, 

8=1 

N-l 

one gets £n = NX — Yl r]^. After this reduction the Poisson brackets between £j, rjj, i,j = 

i=l 

1...N — 1 remain canonical. 
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3.2 Specification of the model 

One may also perform the reduction by the coadjoint action of GLjy acting on the yjl Gaudin 
(ET51) as: 

L Gaudin _^ gL G^ g -\ . A c _^ g A c g- 1 , Vc, g € GL N . (3.11) 
It gives the first class constraint 

n 

^A C = Q (3.12) 

c=l 

and should be supplemented by some gauge fixation \- The reduced phase space is obtained 
by the Poisson reduction 

m Gaudin = 0i x x n //AdG Ljv • (3.13) 

Since dimAdcLjv = N 2 — 1, one gets 

n 

d\mWl Gaudin = ^dim0 c - 2(7V 2 - 1). (3.14) 

c=l 

The SW differential is defined as 

d £Gaudm = ~ dz ( 315 ) 

This reduced model is of our main interest in this paper. Let us start with the example described 
in P. 

Gaudin Model on CP x \{0, 1, q, 00} 

Consider the case of four marked points 0, 1, q, ocd and let A and A 00 be generic orbits of the 
maximal dimension (|3.6p . (|3.7p . while A 1 and A q are those of the minimal dimension ()3.8j) - (|3.9[) . 

4 = ^)> = *y> i<i,j'<JV-i. 



7V-1 

^ = 1> & = E ^i, c^trA 1 

and 

4 = efrj, {e?,»7j} = <%, 1 < *, j < jv - 1, 

iV-1 
i=l 



(3.16) 



(3.17) 



Reduction by AdcLjv leads to (|3.12|) : 

q = A + A 1 + A q + A°° = (3.18) 

with the gauge fixed. The reduction procedure can be done in two steps which deal accordingly 
with the non-diagonal and diagonal parts of the moment map f|3. 18|) . 

The first part of the reduction: 

non-diag(£>) = non-diag (A + A 1 + A q + ,4°°) = , , . 

X = non-diag (^°°) = . K ■ ) 



3 SL 2 acts on CP 1 by rational transformations, which allows one to fix any three of the marked point 
to be 0, 1, 00. 
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i.e. 

A°° = T = diag(ui,...,ujv). (3.20) 

The second part of the reduction involves the rest of the gauge group which is Stab(T) ~ Sj, 
the Cartan subgroup of GLjy: 

diag(e) = diag (A + A 1 + A« + A°°) = , 2 
Xfi = , 

where Xfj are some fixing the Ad^ action. It should be also mentioned here that the spectrum 
of ^4° is fixed 

A^diagOii,...,^), (3.22) 

and therefore 

N 



det( K + T + A 1 +A q ) = m). (3.23) 

Let us calculate the dimension of the reduced phase space. From (|3.14p one has: 
dim m Gaudin = dim O + dim O x + dim O q + dim - 2(N 2 - 1) 

(3.24) 



= 2 dim O max + 2 dim O min - 2{N 2 - 1) 2 N -2 

After the first step (|3.20p of the reduction (|3.13p one obtains a 2 x 2(iV — l)-dimensional phase 
space. Then, the second step of the reduction (by Adg^^/^ ~ Adg, dim Adf D = N — 1) leads 
to the dimension 2{N - 1) as in (g^j]) . 

Below (in Section 5.3) we present the Poisson map of the Gaudin phase space to the Heisen- 
berg chain phase space. We will perform only the first step of the reduction and this gives us 
the exact coincidence of the Poisson structures under the change of variables. 

Reflection symmetry. One can easily see that our Gaudin model (|3.16[) - (|3.22|) possesses the 
following Z2 symmetry: 

(q,z,A°,A\A",T) .— ► z~ l , T, A x ,A q , A ) (3.25) 

Indeed, the transformations (|3.25p do not change the Lax matrix 1-from: 

L G (z)dz = - (T - £ - g^L) & ^ L G {z)dz (3.26) 



or 



L G (z) W -z 2 L G (z) (3-27) 



The described symmetry also remains unchanged in the quantum case since the quantization 
(|6.ip y 1 y d% is in agreement with ()3.27|) : 

d z -L G {z) -z 2 [d z - L G (z)) . (3-28) 

In the gl 2 case, the reflection symmetry structure of this model written in the elliptic 
parametrization [611 [62] was observed in 
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Gaudin Model on C¥ 1 \{z 1 , . . . , z n } 

In the general case, the construction is similar to the previous example. Let z\ = and z n = oo. 
The specific configuration of the Gaudin model under consideration is 

Spec(A°) = (fj,x,...,fi N ) , 

Spec(A 0O ) = (v 1 ,...,v N ), (3.29) 
A U = W • 



(3.30) 



The reduction constraints are 

non-diag(£-) =0, g = A + A 1 + A q + A°° , 
X = non-diag (A°°) = 

at the first step and 

diag(p) = 0, xsj = 
at the second one. The dimension of the reduced phase space is equal to 

dim Wlf audin = (2N - 2) (n - 2) (3.31) 

after the first step of the reduction and finally 

dim m Gaudin = (2N - 2) (n - 3) . (3.32) 

To compute the number of parameters notice that whenever the number of the marked points 
is increased by one this adds two constants (the coordinate of the point and the nontrivial 
eigenvalue of the minimal coadjoint orbit). Then, taking into account ()3.53|) for the space of 
parameters JH, one obtains 

dim (^paudin^ = 2N + 1 + 2(n - A) = 2(N + n) - 7 . (3.33) 

The Lax matrix 

n-l 

L Gaudin (z) = -.4° + V A c 

under the constraint g = can be written in the form 

£Gaudin( z ) = i ^_ T + n^ 1 -^-A c ^j (3.34) 

n-l 

since A = -A°° - ^ A c , A°° = T. 

c=2 

3.3 Poisson reduction 

The Dirac procedure [33J allows one to calculate the reduced Poisson structure in terms of the 
initial brackets on-shell. Let us perform the reduction procedure. There is no any distinguish- 
able way to fix the action \?) (13.21H . This is why we make only the first step (13.301) . 

Suppose we deal with the constraints h = (hi, h m ) = and the matrix of Poisson brackets 
between the constraints on-shell 

c« = {fciA-}L=o ( 3 - 35 ) 
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is non-degenerate at the generic point of the phase space (the second class constrains by Dirac 
[33]). Then, the reduced Poisson structure is given by the Dirac formula. For a pair of functions 
/ and g 



{f,g}o = [{f,g} hi}C^{hj,g} 

In our case, there are 2(N 2 — N) constraints 

h=( gij ,A%), i^j = l...N. 
The matrix C (|3.35[) is of the form: 

a + j3 /3 



C 



n—l n— 1 

, , ,, a = {EA c ,EA c }, (3 = {A°°, A°°} 

P P J c=l c=l 



Therefore, the inverse matrix is equal to 

c- 1 = 



a 1 —a 1 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



and it is not degenerate due to the arguments given in |34j . A direct evaluation leads to the 
following results: 



Proposition 1 For the generic Gaudin model \3. 1\) the reduction corresponding to the first 
step 13.19\) gives the following reduced (Dirac) brackets 



{At p A\ l } D = 8 a Mlh j -AlA l ) 



Qpp ~ Qk k — App + A 1 ^ . 



kp pj 



m-epp+A^p-Al + gu-g^-AZ+A™ -r ea-Qhk-A^+A^ 



^Al^l-Sji) 



+ 



A« kj A%(l-6 ik ) 



(3.40) 



for a,b 7^ n. 



Proposition 2 For the special reduced Gaudin model \3. 29}) the reduction corresponding to the 
first step 113. 19\) gives the following reduced (Dirac) brackets 

{ff, V b j}n = -By (** + E V p-Stp +e ) j , 

6l _ etfii-Su) (3.41) 



{£,i>£.j}D v ._ v ._ gj . +eii 



Vi—Vj-gu+ejj 
for a,b 7^ n. 



The formulae for the generic model (13.40P are also valid for the special reduced model (|3.29p . 
Moreover, ()3.40p follows from ()3.4ip via the initial parametrization A%j = ££rjj. 

3.4 Spectral curve 

Let us again start from the example considered in pp. 
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Proposition 3 [TJ, I13j The spectral curve for the Gaudin model defined by the data &3.16\) - 



113. 20\) has the following form: 

Y Gaud™ z ) : (i + q ) det(y+T) ( 1 + j^r, 1 (y + T)- x e + ^rfliy+rr 1 ? 

(3.42) 

= zP. v {y)+ z^qP^y), 
N N 

Pp(v)=H(v-i*i), P-v=H(y + Vi). 

1=1 1=1 

Proof. 

In order to compute the spectral curve (|3.2p we need the following simple Lemma: 

For any given invertible matrix G € Mat{N) and a pair of N-dimensional vectors £ and n: 

1. det(G + e?? T ) = (1 + ?7 T G- 1 0detG, (3.43) 

2. (G + = G- 1 - l+ri r G - H G- l ^ T G-\ (3.44) 

Substituting ([3^20]) into the Lax matrix flU]) with A c defined by (l3TT6]) - (^23l) . one gets: 

i / 41 Aq \ 

L Gaudin (z) = — - ( T q . (3.45) 

z \ z — 1 z — qj 

Therefore, the spectral curve equation det(y — L G (z)) = can be written in the form: 

det(yz + T — ^x?? 1 — £ 9 x ri q ) = (3.46) 

z — 1 z — q 

Introduce 

y = yz. (3.47) 

Applying (|3.43p . (|3.44p from Lemma 1 twice, one obtains: 

T™(y, z) : det(y + T) (l - ^- lV \y + T)-^ 1 - ^- q r?(y + T)- 1 ^ 

(3.48) 

+ {z -i) (z - q) (riHy+r^e ■ ^(y+V^-riHy + V 1 ? ■ ^(y+T)-^ 1 )) =0 

The expression in the l.h.s. of (|3.48j) contains only simple poles at z = 1 and z = q. The poles at 
y = — vi are apparent. Indeed, it is easy to check that the second order poles are cancelled out 
in the second line of ()3.48p while the factor det(y + T) cancell the first order poles. Moreover, 
let us compute det(y — A ) for ()3.20p - ()3.23p in the same way: 



det(y - A°) = det(y + T)(l + ^{y + T)-^ 1 + rf{y + T)- 1 ^ 

(3.49) 

N N K ' 

+r] ^y + T)- i e ■ ^(y+rr^-vHy + T)- X e • ^(y+T)- 1 ? 1 = n (y - m) 



i=l 



Then, the final answer is achieved by plugging the second line of (|3.49p into the second line of 
(BUS) : 

z 2 det(y+T)-zdet(2/+T) (^(y+T)- 1 ^ 1 + qr] q {y+T)- x ^+q+l) 

(3.50) 

+qdet{y-A°) = 0, 
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N N 

det(y-A°)=H(y - det(y+T) = fj(y + Vi ). 
i=l i=l 

Notice, that the case of the sIn Gaudin model differs from the gljy case by the shift A° 



A c — jrtxA c . Therefore, the spectral curves differ by 



trA 

y ->y 



1 trA c 

+ w E • (3 ' 51) 

%e{o,i,?} 



The space of parameters of the Gaudin spectral curve is described by the following set: 

^Gaudin = {vu ^ VNtfli> ^ trA l 5 trA9j g} _ (3 52) 

Taking into account the possible shift of y, the number of independent parameters is equal to 

dimm Gaudin =2N + 1. (3.53) 
The case of arbitrary number of marked points (|3.29p is considered below (see Theorem 2). 

3.5 Simplest example: gl 2 on CP 1 \{0, 1, q, oo} 

Let us now calculate the spectral curve. In the gl 2 case 

det ^2 A i = ^2 det Ai + ^2 trA i tlA i ~ ^AiAj, A* G gl 2 . 

i i i<j 

The special case (|3.16p - (|3.22p for gl 2 means that one deals with four 2x2 matrices 

A ~ diag( Ml! /x 2 ) ; A 1 = e x r, 1 , A« = x v *, 

where £ , £ 9 , rj 1 , rfi are 2-dimensional vectors and A 00 = diag(i>i,U2) with the condition (|3.20p . 
For the spectral curve we have 

(3.54) 

+ pr ^det I + q — {z _ 1)(z _ q) Q — q — ) - U, 

where in the second line we used A = -T - A 1 - A q . Alternatively, one can make the shift 
(|3.5ip y —> y + ^ (j§4f + Q^~^ + m ^ 2 ^ , which corresponds to the traceless case gl 2 — > sl 2 , i.e. 
Ai ~ diag(^j, —Vi). In this case the spectral curve can be written in the following form: 

f - Viz) = -, -H (3.55) 

where H is the Hamiltonian function on the phase space 

// ,J< ( • q * M - f 2 " 

and the potential reads 



^(^1+^2 + ~~\ UX I ( 3-56 ) 



K (Zj " + (js _ i)2 + {z _ q) 2 z{z _ 1} i^- 57 ) 

with v = vi = ±tiA\ v q = ±trA*, i/oo = *a^a 

Remark: The coset space ()3.13j) in this case is the phase space of the Painleve VI equation 
[65] in the Schlesinger description [66] . 
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3.6 AHH duality 



In [29j M.R.Adams, J.Harnad and J.Hurtubise suggested a duality between the classical Gaudin- 
Schlesinger models of different ranks and numbers of the marked points. Their description of 
the Gaudin model differs from ours by the constant term Y in the Lax matrix (or connection 
in the isomonodromic case): 

M A c 

L G ahh{z) = Y + Y.^—' Y = di ag ( yi ,...,y N ), A c eg\ N . (3.58) 

c=l 

The difference is essential, since 7/0 leads to appearance of the second order pole at oo 
for L^ HH (z)dz. The phase space is also different in this case. It is a direct product of the 
coadjoint orbits (equipped with a natural Poisson-Lie structure) factorized by the stabilizer of 
Y: (O 1 x ... x O m ) //Stab(Y). 

In the case when all A c are of rank 1, the dual Lax matrix is the gl A /- valued function with 
Y = diag (zi, ...,zm) and TV" marked points at y\, ...,i/n'- 

N A c 

L < Ahh(z) = Y + T , Y = diag(z l ,...,z M ), A c G gl M . (3.59) 

The duality implies the relation between the spectral curves: 

det(Y - z) det(L% HH (z) - A) = det(Y - A) det{i% H {\) - z) . (3.60) 

Dimensions of the phase spaces of the both models are equal to 2(iV — 1)(M — 1) and the 
number of parameters is 2{N + M) — 3. Indeed, for (|3,58p dimension of the unreduced model 

M 

is dim^4 c = M x 2(A^ — 1). The reduction by the coadjoint action of the Cartan subgroup 

c=l 

Stab(y) ~ Sj C GL N , dim Ad fi = N - 1 leads to 

M 

J^dimA c - 2dim(Stab(y)) = M(2N - 2) - 2(iV - 1) = 2(N - 1)(M - 1) . 

c=l 

4 Heisenberg Chain 

4.1 GL2 XXX Heisenberg chain 

Let x be a local coordinate on CP 1 . Define the Lax operators as a set of GL2-valued functions 

Li{x)=x-Xi + S\ i = l...N (4.1) 

3 

where S l = ^ S % a a a are matrices from skj and {x{} is a set of points on CP 1 . Each Lj(x) is 

a=l 

assigned to the i-th site of one-dimensional lattice. Then the transfer-matrix is defined as 

T{x) = V(q)L N (x)...L l (x) (4.2) 

where the "twist" V(q) is a constant GL2-valued matrix. Following [46] we choose V depending 
on a parameter q 

' 1 



V = I W I , q = const (4.3) 



17 



with eigenvalues and j^. The spectral curve and the SW differential are 

f Hoiscn (u>, x) : det(w - T[x)) = 0, (4.4) 

dS xxx = x dw 
w 

Expanding the determinant of the 2x2 matrix, one obtains 

w -trT(x) +W' 1 detT(x) = 0. (4.6) 
The Poisson structure consists of N copies of the Lie-Poisson sl2-brackets at each site: 

{S k a ,S l p } = S kl V^ls a ^S k . (4.7) 
These brackets are generated by the quadratic r-matrix structure: 

{L(z) ® L ( w )} = [ r ( z , w ),L(z) ® L{w)\ , 

'•(«.«') = ijfl®l+EiT a ®ffaj- ^ 
The values of the Casimir functions are defined by the eigenvalues Spec(£ l ) = (Ki, —Ki): 

3 

Kf = J2 S aSi, i = l...N (4.9) 

a=l 

Thus, there is only a pair of independent variables at each site. Since det Li(x) = (x — xi) 2 — Kf 
the spectral curve (j4.6j) is now written in the following form 

w - tvT(x) + w' 1 . q xo K + {x)K-{x) = 0, (4.10) 
(1 + q) 1 

N 

K ± (x) = \\(x-mf) mf = x i ±K i . (4.11) 
i=l 

At this stage the phase space of the model is 2A r -dimensional 

9Jt Heison = {S l a , i = 1. . .N, a = 1. . .3 | K { = consti e C} , dim«m Heisen = 2N (4.12) 
Notice that the functions K (x) are constant on J0t Hcisen ) while 

N 

trT(x) = x N + ^' l Hf ciscn , (4.13) 

where |iJP eisen } i s a set of N polynomial functions of degrees N+l — i on 5Q? Heisen which mutually 
Poisson commute {ilf eisen , Hf ei3en } = with respect to (|47f|) due to 

Let us also get in touch here with the parameters of the four dimensional gauge theory: 
mf correspond to the masses of the 2N hypermultiplets in the fundamental representation, 
q = e 27 " T i s related to the ultraviolet value of complex coupling constant r and xi parameterize 
the vacuum expectation values of the scalar fields. 
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4.2 Reduced phase space and spectral curve 



The Hamiltonian Hff C1SEn is linear: 



rrHcisen 

N 



/ N \ N 
V i=l J i=l 



It means conservation of the projection of the total "spin" along the "vector part" of V: V = 
^tr (Va). Obviously, it is generated by the adjoint action of g G GL2 which keeps V unchanged, 
i.e. g £ Stab(V): 

g € GL2 : S'^gS'g' 1 , V^ g Vg- l = V (4.15) 

By analogy with the usage of the center of mass frame in many-body systems (when the 
total momentum is conserved) we are going to resolve the equatiorQ 

^Heiscn _ ^jgj 

with respect to any variable, say , and to reduce the phase space 9Jt Hciscn by 2 dimensions: the 
constraint (|4.16|) should be accompanied by some "gauge fixing" constraint, which corresponds 
to a choice of nontrivial action (|4.15p on S . Let us denote this second constraint as Xn{S n ) = 
0. 

It is easy to show that the Poisson brackets (j4.7|) between S k , k = 1...N — 1 are not changed 
by this reduction. Finally, the reduced phase space 

m HciBcn = {S^ i = 1...N, a = 1...3 I Ki = const;, F^ ciscn = const, Xn{S N ) = 0} 

(4.17) 

dim2Jt Heison = 2N - 2 

has dimension 2N — 2 and the integrable dynamics is generated by Hf clBon ,...,H^ c ^ n with 
S N = S N (S 1 ,...,S N ~ 1 ) on-shell. 

Let us also rewrite the spectral curve (14. 4| . Substituting w = w j^K + (x) into gSD, one 
gets 

r Hciscn (u;, x) : trT(x) - -^-wK + {x) - w~ l K-{x) = (4.18) 

or 



x 



N + ^A^=l-^h) [w 2 AK+(x)-w(l + q)AtiT(x)+qAK-(x)]=0, (4.19) 



where A/(x) = f(x) — x N . The ultimate curve depends on N — 1 Hamiltonians and on the 
following set of 2N + 2 parameters: 

^Heiscn = ^ ^ R ^ ^ ^ tf^ g } ( m ± m ± , f/^ Q } . (4.20) 

In fact, one of {xi} can be moved away by the shift of x. Therefore, the number of independent 
parameters is equal to 

dim$n Heison = 2N + 1. (4.21) 



4 The natural choice from the point of view of the gauge theory is to set = because this eliminates 
the U(l) factor and leaves the SU(N) theory. However, we keep the constant arbitrary. 
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4.3 Simplest example: 2-site chain 

Let 

trT(x) = x 2 + H 2 x + H 1 , K± = (x - mf )(x - mf ) (4.22) 
Substituting it into (I4.18j) with the change of variables x = yz, w = z, one has 

y 2 z 2 + H 2 yz + Hi = —^—(yz - mf)(yz - m%) + (yz - m^)(yz - mj) (4.23) 

This equation is easily reduced to the form 

f - b(z)y + c(z) = 0, (4.24) 

with 



which is rewritten in the form 



y 2 - V(z) = , l — £ -H (4.27) 

V ' z(z-l)(z-q) 



where H is a Hamiltonian (linear in Hi) and 



with 

m = 



TOq = —L-x — -, TO 



(4.25) 



t / s _ (q+l)H2+q(m+ )+(m 1 +m 2 ) (m 1 +m 2 ) (g+l)JJ2+(mj"+m 2 ^)+g(m 1 +to 2 ) 
°\ Z >- _{q-l)(z-q)_ + z (q-l)(z-l) 

r \ g+1 m i m 2 ~Hi m i m 2 +<Z m i~ m 2~ — (l+?)-Hi q m i m 2 +'m-~^m J — (l+q)Hi m 1 rn 2 

C y Z > - — z + q(q-l)(z-q) (g-l)(^-l) + ' 

Making a shift y — > y + \b(z), one comes to y 2 + c(z) — \b 2 (z) = 

y 2 + c(z)-^b 2 (z) = (4.26) 



= £ + ^ + ^ - ^jffi"'- (4.28) 



2 

m i = 27^TT ((9 + 1)^2 + to+ + to+ + q(m^ + rrq)) , (4.29) 



m Q = ~2(F+) ((? + 1 )^2 + to x + to 2 + q{rn^ + mj)) . 

This should be compared with (|3.55|) - (|3.5Tj) . Note that H 2 is a constant (|4,16p . Therefore, mi 
and to 9 are constant as well since they are independent of Hi. The choice of signs in (I4.29P 
does not follow from (|4.28p . Our choice of signs will be justified in Section [6l 



4.4 Higher rank chain 

Let us now consider the GL^ model 

L i (x) = x-x l + S i , ^Gghfc (4.30) 

with the transfer-matrix 

T(x) = V(q)Li(x)L 2 (x)...L N (x), (4.31) 

where V(q) = diag(i?i, ...,Vk) is a diagonal matrix (its entries will be determined later) and the 
Poisson brackets are 

{L n (x) ® L m (y)} = 5 nm [r(x, y),L m (x) <g> L m (y)] . (4.32) 
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The r-matrix is defined as 

k 

r(x, y) = , Pi 2 = V £ a6 ® P ba . (4.33) 

cc — W 

w a,fe=l 

Then 

^} = (S'aJcb - Si b 5 ad ) . (4.34) 
By virtue of the r-matrix relations (|4.32p . one has the same brackets for the transfer-matrix: 

{T(x) ® T(y)} = [r(x, y),T(x) T(y)\ . (4.35) 

Therefore, 



{Tab{z),T cd (w)} = (T cb (z)T ad (w) - T ad (z)T cb (w)) , 

{T ab (z),T cd (z)} = T ad (z)d z T cb (z) - T cb (z)d z T ad {z) . 
Let us now represent the transfer-matrix in the form of a sum: 

\ \ v , x -v / N & 



(4.36) 



T(x) = \{(x - x 3 ) V{q) [J ( 1 + — J = 11^ - *i) +E^7 ■ ( 437 ) 

j=l i=l l ' j=l \ i=l / 

The gl fe -valued coefficients S l are the residues of the expression T{x)/\\{x — Xj), i.e. 

j 

~ s i n*i) _ 1 / T(x) dx 

Uj^i(xi - xj) 2m J y . ll ( (.r - arj) ' 

where 7« is a small contour encircling x%. Let us now calculate the brackets for the new variables 
S l in the spin chain. 



Proposition 4 The Poisson brackets for S l ^4-38 ) are of the following form: 



qi qj _ qi qj I qi qk _ qi qk 

S'J = (1 - ch ad _ ad A + & S ad V cb - S l cb V ad + £ ad Cb _ Cb ad I • (4-39) 



X{ Xj y ^ ^ 



Proof. 



For i 7^ j the brackets are easily obtained from the r-matrix relations for T(x) (14.361) : 

{S ab , S J cd } = ~ S *d S tb for • ^ ■ (44Q) 

For i = j the arguments of the transfer-matrices coincide. One can use either local expansion 
of (|4.36p or integral representation (|4.38p to overcome this problem: 



{S ab , SU) = L L dx dy{T ab (x),T cd (y)} U ^ 



l l 



4 f ^ sffi^a p ^-^ = t„ * ™ Srn^f • ■ (441) 
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Now we shift the contour 7^ to infinity so that it transforms to small contours encircling Xj 
for j ^ i and the contour around x = 00, all of them going counterclockwise. The integral is 
computed by residues: 

r qi qi 1 _ X r j^ T ^b{ x )S ad -T ad {x)S l cb 

i\b^cdi / ; A >.r r,.||..r r., 

_ X ^^ T cb(x)S ad —T ad (x)S l cb _ ^ ScbSj d -Sg d S l cb _ -tr ni _ t/ ,01 _ 
= ^ad^cb — S l cb V a d + ^ j. jj, (Sad&cb ~ ^cb^ad) ' 

Though the formulae given above are valid for generic S l , we consider the case when the site- 
variables are the coadjoint orbits (of GL&) of the minimal dimension (see (|3.8p ) in this paper: 

dim(^) = 2(k - 1) , Vi = 1...N . (4.43) 

Then the dimension of the phase space is equal to 

dim (9JT Hcison ) = N(2k - 2) - (2k - 2) = 2{k - l)(N - 1) (4.44) 

where " — (2k — 2)" comes from the reduction by the Cartan subgroup of GL& (14.171) . To 
compute the number of parameters note that, with increasing the rank by unit, one adds two 
more constants (the component of V and the value of the Hamiltonian corresponding to the 
action of the Cartan subgroup of GL/u). Then, taking into account (|4.2ip . one obtains 

dim (m Hciscn ) = 2N + 1 + 2(k - 2) = 2(N + k) - 3 (4.45) 

for the dimension of the space of parameters. 



5 Classical Duality 



Let us first recall the result of pp. 



5.1 Duality for GL2 chain 

As it follows from Proposition 3 (Section 3.4), the spectral curves of the special reduced Gaudin 
model (I3.16p - (13.22p and the XXX chain have the same forms ()4.18|) and (|3.42[) . Moreover, the 
dimensions of the phase spaces (|4.17j) . (|3.24p and the spaces of parameters (|4.20p - (|4.2ip . (|3.52p ~ 
(|3.53p can be identified as well. This is summarized in the following 



Theorem 1 The N-site GL2 Heisenberg XXX chain defined by Iji4-1\ )- [4~£1\ ) an d the gl N Gaudin 
model \3.11\) - f3.53\) are spectrally dual at the classical level with the change of variables: 



I y 



w, 



x_ 

w 



(5-1) 



the following identification of the parameters ^4-20 ) and \3.52\) : 



Tu- 



rn; 



1 < i < N, 



1 < i < N, 



(5.2) 



N 



k=l 
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and the following relation between the generating functions of the Hamiltonians: 

trT^(y) = det(y+T) (l+ 1 ^(y+T)-^ 1 + ^(y+T)" 1 ^) (5.3) 



The simplest example of the duality comes from comparison of (I4.27p - (|4.29p with (13.550 - 
(1337]). 

5.2 Duality for higher rank chains 

All the statements of the previous sections work not only for the GL2 spin chain and the 4-point 
Gaudin system, but also for the GL^ spin chain and the k + 2-point Gaudin model. Indeed, 
from (|3.32p , (|3.33|) and (|4.44|) , (|4.45|) one concludes that the dimensions of the phase spaces and 
spaces of parameters are equal for the Gaudin model (|3.29j) - ()3.33[) and the Heisenberg chain 
([Q4>([Q4l if n = k + 2. 

Theorem 2 The N-site GL& Heisenberg XXX chain ( |^.37| )- |^.^5[ ) and the gl/y Gaudin model 
on CP 1 \{zi, . . . ,z n } h3.29\) - F3.33\) are spectrally dual at the classical level with 

1. the change of the spectral problem variables: 

z = w , 



(5.4) 



v = — 

& id 



2. the following identification of the parameters: 

n = k + 2, (5.5) 

V = diag(z 2 , ...,z n -i) G gift , (5-6) 

A99 = r ti = Vi = - Xl . (5.7) 

3. and the following change of dynamical variables: 
V-W) =&r$, a, 6 = 1.. .k, i = l...N, (5.8) 

/ ab 

Proof. 

Let us see explicitly how the identification between the two systems manifests itself at the 
level of spectral curves. One writes for the spectral curve of the spin chain: 

de,Ln(»- Ij y(l + |M.)to. (5.9) 



For the k + 2 point Gaudin system: 



det (yz - L G ) = det ( yz + - V Z -^L1 \ = . (5.10) 



NxN NxN \ z — ' Z — Z„ 
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We denote yz as x and assume as in the preceding sections that the matrix is diagonal with 



the diagonal entries v% 



1 . . . N. Then one has: 
detTVxTV (x + Aoo-J^ 



k 

a=l 



N 



f] (x + Ai) det NxN 
i=l 



Z-Za 



i-(, + ^)-x^ =l -5gE 



(5.11) 



Let us introduce D = (x + A^) 1 Ea=i 



k z a ( a (r 1 a )' 1 



Expanding the determinant into powers of 



-D, one gets the sum of terms, each term being a product of traces of D n . A typical contribution 
looks like c nij ... jnj trD ni ■ ■ ■ D n K Let us point out that the coefficients c nii „. jVlj do not depend 
on N, the size of the matrix Lq- We are going to prove therefore that each trace of D n can 
be rewritten as the trace of the n-th power of a k x k matrix. This will prove the equivalence 
between the two spectral curves. 

The derivation goes as follows. Consider 



^NxnD 

N 

= E 



1 



N k 

E E x+Ai ' {z-z ai ) 

il...i M =l ai...a,M=t 



7 f il n 2 
Sa^ 'lai 



X+Ai 



(Z-Z a ,.) 



1^ tr kxk [{z-Z) Z {x+A ^ ■■■{z- Z) Z {x+Aim) 



(5.12) 



t r fcxfc 



■ N 

E 



(z-z)-^z V '(ey 

x+Ai 



M 



where the traces in the second and third lines are taken over the indices aj and Z is the 
gl fc -valued diagonal matrix with the diagonal entries z a : 

Z = diag(zi, ...,z k ) . 



Thus, one has: 
detjVxAr [x+A 



Ek 
a=l 



za^ a (v a y 

Z-Z a 



N 

J] (x+Ai) det fcxfc 
i=i 



N {z~Z)- 1 Zr 1 i (e) T 
x+A t 



n (p+Aj 

fei det 

II {Z-Za) 



kxk 



N r, l {c) T 

x+Ai 



(5.13) 



The spectral curves coincide provided one identifies Z = V[j, S l = V^(r] l ) T and V{ = — Xj. ■ 

Remark: in fact, any Gaudin system can be rewritten as a spin chain by merging together 
several points containing minimal rank orbits. This will correspond to merging of the eigenvalues 
of the matrix V in the spin chain. Conversely, any spin chain can be rewritten as a Gaudin 
system by analogously merging several of the points Xj. The matrix A°° in the Gaudin system 
then has several coinciding eigenvalues. 



5.3 Poisson map 

In this paragraph we show that the identification (|5.8j) provides a Poisson map between the 
two models. Our strategy is the following: we start from the Gaudin model (|3.29j) - (j3.34p with 

5 From now on, we use the different overall normalization for z and, hence, for V: that is, for k = 2 
we set V — diag ( — j— , J instead of V = diag (1, q). 
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the quadratic Poisson brackets (|3,4ip . Then we apply the AHH duality transformation (|3.59p 
to 2 L Gaudm . After that we show that the gl n _ 2 -valued residues of the dual model obey the same 
Poisson structure as S l f|4.39j) . 



Denote 



^ 6 = «egl fc . (5.14) 

In fact, the proof of Theorem 2 follows by applying the AHH duality (|3.60p . One can define 
the AHH dual of the Gaudin model (ET34"j) 



* din (z) = - I -T + V -^-A c , T = diagK v N ) 

Z\ ^*"*c J 

following the recipe ()3. 59 j) : 

r -i N 

ARR\zL G ^(z)\(x) = Z+ £ £JL (5.15) 

^ ' i=l ' 

However, the Poisson structures of both models are quadratic in our case (instead of the linear 
brackets in AHH [29J). 

Theorem 3 Under the AHH duality transformation the phase spaces of the gl N reduced k + 2- 
point Gaudin model \3.29\) . 13.31\) . j3.41\ ) and the XXX Heisenberg chain ^4.3^ - [f3^ are 
related by 115. 8\) 

(V 1 !?) =Bl b , a,b = l...k, i = l...N. 

This map is Poisson, i.e. the Poisson brackets for B l ^5.1$ generated by the Dirac brackets 
{3.4 1\ ) are of the same form as for V^S 1 /[4-39\ ): 

{BL,Bl d } = (1 - 8^) Vi _ v .} eii+e .. {Bi d B^ - B* cb B ad )+ 

(5.16) 

R* A, _ R* A , _|_ V , Q-d a cb D cb a ad j 

with 

Xi = Qu - Vi . (5-17) 
The variables qu and Vi are the Casimir functions in the AHH dual Gaudin model: 

{v i ,A b kl } D = {v l ,A b kl } = forb=l...k, 

{eu,A b kl } D = { gii , A b kl } = for b = 1 . . . k . &) 

The proof follows directly from the Dirac brackets (|3.4ip . Notice that qu are the Casimir 
functions. Therefore, one can put qu = and reproduce (|5.7p . Then, comparing (|5.15p and 
(|4.37p . one gets the exact equivalence of the models: 



AHH|zL Gaudin (z)|(x) = T Hciscn (j;)n • 



(5.19) 
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6 Quantum Duality 



6.1 Simplest example 

Let us start again with the simplest case (sl2 4-points reduced Gaudin model and 2-cite chain) 
corresponding to the four-dimensional Yang-Mills theory with SU(2) gauge group. The Baxter 
equation for the XXX spin chain is known to lead to the Bohr-Sommerfeld integrals that describe 
the SW system giving rise to the right Nekrasov prepotential [TT]. This Baxter equation can 
be thought of as a quantization of the corresponding spectral curve ([4.230 : 



trT(x) = yV + H 2 yz + H 1 



1 + 



1 K+(yz)+- 1 



zl + q 



K~{yz), 



where 



K + (x) = (x — mf)(x — mjj"), K (x) = (x — m 1 )(x — m 2 ). 
The "quantization" implies that (for x = zy) 

y — > hd z 



(6.1) 



At this stage we need to fix the ordering. We use the normal ordering putting d z to the right 
of all z's. This provides the following rules 



K+(x) :=: K+{x) : +hx, 
jf-(x) :=: K~{x) : +hx, 
trTji(x) :=: trT(x) : +Hx. 



This recipe follows from 



x 2 Qaunt R 2 Z Q Qz = h 2 z 2Q2 + h 2 z Q z = .j,2. +h::i , : 



After substituting ([6.10 . one has 

[h 2 dl - b h (z)hd z + c h (z)] ipHc iS cn(z) = 0, 

with 

b h {z) = b(z) - §, c h (z) = c(z), 
where b(z) and c(z) are given by ([4.250 : 



(6.2) 

(6.3) 

(6.4) 
(6.5) 



b(z) 
c(z) 



(q+l)H2+q(mf )+(m 1 +m 2 ) (m 1 +m 2 ) (q+l)H2+(m+ +m^)+q(m 1 +m 2 ) 



(q-l)(z-q) 



+ 



(q-l)(z-l) 



q+1 m 1 m 2 —Hi m 1 m 2 +qm^m 2 — (l+q)H\ qm 1 rn 2 +m~^m 2 — (l+q)Hi m 1 m 



q : ' q(q-l)(z-q) (g-l)(z-l) 

The "classical" shift (j3.5ip y — > y + \b(z) is replaced here by the change 



+ 



'0Heisen( 2 O 



V ; Gaudm(^)exp I i [dz b h (z) 



z 



which leads to 

h 2 d 2 z + c n (z) - -b\{z) + -(d z b h (z)) Veaudmv 
Making the same calculations as in the classical case, one finally gets 



0. 



h 2 d 2 z + v H (z) 



z(z - l)(z - q) 



Hn 



IpHeiseniz) — 0, 



(6.6) 



(6.7) 



(6.8) 
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where H% = H + hm q and 



M(„\ — m \~h: i mi (mi -ft) | m q (m q ~h) ml+m±(m±-h)+m q (m, q -K)-m 2 oa (a q\ 

Vh[Z) - Z 2 + (z _l)2 + {z _ q) 2 ^±15 

mo, mi, m g and m^ being given by fl4.29j) . Equation (|6.8p is exactly the one derived from 
consideration of the conformal blockqj. In [12J it was shown that this equation gives the 
quantum periods (monodromies around the A and B cycles) from which one indeed obtains the 
correct Nekrasov functions. The two wave functions ^ He isen and V'Gaudm give slightly different 
periods and thus different prepotential because of the factor bn in (|6.6p . The difference for the 
period integrals is 

An = <b b h (z) dz . (6.10) 



However, it does not affect the physics for the following reasons. Indeed, as we briefly mentioned 
in Section I4.2l the linear Hamiltonian Hi = ^2 k which enters bfr is set to vanish for the SU(iV) 
gauge theory. The remaining terms do not depend on the Hamiltonians Hi or equivalently on 
which in the gauge theory correspond to the vacuum moduli. Thus, the contribution All of 
b% to the monodromies is a constant (depending only on the masses nit), its contribution to the 
prepotential is linear in vacuum moduli and does not enter the low-energy effective lagrangian 
C = J d^9F{^) due to the identity Jd 4 9^ = 0. 

Formulae (16. 6ft and (|6.7p present the relation for the sl2 Gaudin model while the gl 2 case is 
much simpler: 

^Hciscn(z) = V'Gaudin(^) • (6-H) 



6.2 General case 

Quantization of the XXX chain spectral curve (|3.50[) with the SW differential (|4.5p means that 
x should be simply replaced by hwd w . For example, in the GL2 case one gets the Baxter 
equation: 

tr T(hwd w ) - -^K + (hwd w ) - - — ^—K-fHwdS) * Hciscn (w) = 0. (6.12) 
l + q (l + q)w J 

Equivalently, for the Gaudin spectral curve (|3.42p the quantization is given by the replacement 

y ->• hd z -. 



n (zHa t +vi) + E ,ik x + ; n + «o- 

\ i=1 k=1 ' ^ k (6.13) 

-^TT Ft (^d z + Vi ) -z-^U (zHd z - ih) Hf^(z) = . 

i=l i=l J 

Obviously, the differential operators in the brackets of (|6.12p and (|6. 13j) can be identified in the 
same way as the classical spectral curves did. 

The Baxter equation (|6.12p looks similar to the classical equation for the spectral curve 
except for the substitution x —> hzd z . All the differentials zd z are placed to the right of all 
the functions of z (which is consistent with our previous normal ordering under the replace 
z — > log z). 



6 Different conventions on what to call the mass parameters in the gauge theory after the e-deformation 
exist. To compare our expressions with [5] one needs to perform the shift: m^ ur = m^ GT ± h/2. 
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To obtain a similar expression for the quantum Gaudin system one also needs to place all 
the differentials zd z to the right. After doing this, the equivalence is evidently analogous to the 
classical case. Thus, the definition of the quantum "determinant" must be the following: 

" det " (hzd z + zL Gaudin (z)) d = all zd z to the right. (6.14) 



Theorem 4 The N-site GL& Heisenberg XXX chain and the gl^r k + 2-point Gaudin model are 
spectrally dual at the quantum level with the following relation between the wave functions: 

IpXXx(z) = ^GaudJz) , (6.15) 

i.e. the Baxter equation of the Heisenberg chain can be rewritten as the quantization of the 
Gaudin model spectral curve with the ordering $[6.14\ )- 

A small problem arises only for the sIat Gaudin model. In this case, one has an additional 
factor, just as in (16.61) . However, this factor is controllable. Let us compute it in the case of 
GL2 chain, i.e. in the case of 4-point Gaudin model duality. 

The Baxter equation for the XXX spin chain can be rewritten as a polynomial in d z (instead 
of zd z ). To do so, one needs to make the ordering in (zd z ) n . To this end, let us represent the 
action of (zd z ) n as a part of the dilatation operator 



, e azd *f(z) = f(e a z). (6.16) 

a=0 



Then, using the Taylor expansion for f(z + z(e az — 1)), one gets 



k=l 1=1 v >' ' 1=1 ■ \ / 

The coefficients have the following properties: 

<5J = 1, ^1 = 1^1, C£ = l. (6.18) 
For instance, the latter one is the well-known combinatorial identity [63 : 



1 — n V / v 



n-k ( n \ ,„r _ / 0, if r < n, 
nl, if r = n. 



k=0 

Finally, for K + (zd z ) one has 

N N n 



K + (hzd z ) = Y,(-l) n Vn({m + })(frzd z ) n = J2(-l) n c- n ({m + })h n C*z k d* (6.19) 



n=l n=l k=l 



and similarly for K , where cr n ({m + }) are elementary symmetric polynomials, and trT. Now 
one can rewrite the Baxter equation in the desired form in which all d z are placed to the right 



52 J2 (Hn+l - ^-^a n ({rn + }) - ^-^a n ({m-})) h n C k n z k d k z rfe = (6.20) 



n=l k=l 
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Analogously to the transition from ()4.18j) to ()4.19j) . one may write the equation in the form: 



iV-2 



h N a? - bnWtF^d?- 1 + c m {z)h m d^ U = o (6.21) 

V m=0 / 

Prom fleSH) , (pUBD and (IBTTHI) it follows that 

iB(2) .^^(^ + £%< + i£|^)_ ft ^). {6 , 2) 

(z - l)(z - g) ^ 1 + q {l + q)z ) 2z 

Again, one may eliminate b^(z) from the equation. By the construction, we come to 

Theorem 5 p] The N-site GL2 Heisenberg XXX chain defined by ^4.1\ )- ^4^l\ ) and the g\ N 
Gaudin model i3.11]) - [3. 53\) are spectrally dual at the quantum level with the following relation 
between the wave functions: 

^xxx(z) = ^ G auMe^ rbhdz . , (6.23) 

i.e. the Baxter equation of the Heisenberg chain can be rewritten as the quantization of the 
Gaudin model spectral curve with the ordering \6.14\ ). 

Indeed, all the classical statements work with our choice of ordering (|6.14p for the Gaudin 
system. Thus, the two systems are exactly equivalent. 



7 Comments and Discussion 

In this section we discuss possible generalizations and relations of the models and the corre- 
sponding duality. First of all, there is a wide class of the monodromy preserving equations 
which can be considered as some generalization of the Gaudin model. 

• Painleve-Schlesinger Equations. In the gl 2 case, the coset space (I3.13j) is the phase 
space of the Painleve VI equation |65| . The Painleve equations can be interpreted as 
equations describing isomonodromic deformations [66] of the linear system 

(d z -L G (z))*(z,q) = 0. (7.1) 

The "deformation" generates dynamics with respect to the marked points. In the case of 
the 4-points Gaudin model, it is only q: 

A q 

(d q -M^(z))^(z,q)=0, M°{z) = — . (7.2) 

z q 

The zero-curvature condition for the above equations is the monodromy preserving equa- 
tion: 

d q L G {z) - d z Mf{z) + [L G (z), M G (z)} = 0. (7.3) 

The residues at 0, 1, q provide the finite-dimensional equations on A c which are called the 
Schlesinger equations: 

d q A° = -l[A°,A% 

d q A^ = ^- q [A\A% (7.4) 
d q A* = \[A°,A*]- ^[A\Ai]. 
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After the reduction f)3. 13j) the later equations reduce to the single one: the Painleve VI 
equation. All the same can be done in the general gl^r case with any number of marked 
points. In [13] Y.Yamada0 considered the higher Painleve equations corresponding to the 
configuration (I3,16l) - (]3.22p . Note that if one replaces d q by dt in all the above equations, 
(|7.4p describes the flow in the unreduced Gaudin model, while the reduced one is just the 
autonomous version of the corresponding Painleve equation(s). The Schlesinger system 
is non- autonomous version of the Gaudin model: 

Mi = %^ -> V^ = %^ (7-5) 

Z% Zj Z{ Zj 

Relations between the Gaudin and Schlesinger systems can be summarized as follows: 

(7.6) 



(A = 






-L{z))W = 


(d ti - Mi{z))V = 




l (d Zl 


-Mi(*))tf = 


d k L=[L,M i ] - 


> d z L 


- d z M -- 


= [L,M t ] 


{H i ,H j } = -> 




- V* 


= {H i: H j } 



(7.7) 

Relation between the Schlesinger system and the Painleve-type equations is given by the 
reduction procedure: 

Unreduced Gaudin — > Schlesinger System 

I I 
Reduced Gaudin — > Painleve Equations 

Let us list open problems arising from the Gaudin-Schlesinger correspondence: 

1. The recipe (17. 6p relates the isomonodromic deformations with the quantization prob- 



lem. The corresponding "isomonodromic quantization" was studied in [68]. The natural 
quantization is given by the Knizhnik-Zamolodchikov-Bernard (KZB) equations for the 
Gaudin model. We hope to find a relation between solutions to the Baxter equation and 
the corresponding KZB equation. 

2. In [62] it was shown that the classical linear problem at the r.h.s. of (|7.6p leads to the 
non-stationary quantum problem for the reduced model written in spectral variables. At 
the same time, V.Bazhanov and V.Mangazeev showed [73] that the Baxter equation of 
some special chain is compatible with the Lame equation. Therefore, one can expect that 
the Baxter equations are compatible with some kind of the non-stationary Schrodinger 
equation. 

3. In [13] the higher (coupled) Painleve-type equations related to the specific configuration 
(|3.16p - (|3.22p were studied. In the simplest case the underlying autonomous system is 
the BC\ Calogero-Inozemtsev model, while its isomonodromy version is the Painleve VI 
equation in the elliptic form [ST]. For the gLy case one can also hope to find an elliptic 
integrable ("spinless") system of Calogero-Gaudin type and the corresponding Painleve- 
type system with N — 1 degrees of freedom and 2N coupling constants. 

4. The elliptic form of the Painleve VI equation can be presented as a non-autonomous 
version of the 1-site XYZ chain interacting with a constant external magnetic fielcH, i.e. 
with dynamical boundary conditions satisfying the reflection equation [SI]. Therefore, 
the spectral duality provides some relation between the XYZ and XXX models which we 
hope to investigate. 



7 This type of equation appeared earlier within the context of the UC-hierarchy in [ST] . 
8 The corresponding mechanical system is the Zhukovsky-Volterra gyrostat. 
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• Quantization. It should be mentioned that we do not impose any boundary conditions 
which provide a valuable quantum problem, i.e. we do not specify wave functions ex- 
plicitly. We compare the Baxter equations which are the quantizations of the spectral 
curves written in separated variables. Alternatively, one can specify the spaces of so- 
lutions initially and then verify their identification through the duality transformation. 
This is the recipe of [69] where the authors considered very close problem in terms of the 
Bethe vectors. Recall, that the Lax matrix of the n-point reduced Gaudin model is of 
the form 

n— 1 / n \ 

L%) = Ei^ = | -T+Eif , (7.8) 

c=l V c=2/ 

where z\ is set to zero and T is the residue at oo. The expression in the brackets (at the 
r.h.s. of (|7.8p ) is similar to the case of Gaudin model considered by E.Mukhin, V.Tarasov 
and A.Varchenko. In [69] they conjectured a duality between the XXX chain and trigono- 
metric Gaudin model at the quantum level, which relates the corresponding Bethe vectors. 
The derivation uses the space of quasi-polynomials (and quasi-exponentials), i.e. the ex- 
plicit form of solutions to the Baxter equations. It is based on the linear (Lie algebra) 
commutation relations, while the Gaudin model under consideration here is the rational 
(although reduced) one. The Poisson structure (which is discussed at the classical level) 
is quadratic. Therefore, one can expect to find a relation between the results of this paper 
and those of [69] by use of a quantum version of the Dirac reduction. At the classical 
level, the main difference comes from the definition of the SW differential. In our case 
dS = Adz with the change of variables to the XXX chain: z = w, A = x/w. In [69J 
the differential (presumably) is dS = A^ and the corresponding change of variables is 
z = w, A = x. Notice also that the expression in the brackets (at the r.h.s. of (|7.8p ) 
contains the constant term (T), which (without the factor 1/z) leads to the second order 
pole at oo for L G (z)dz. This type of Gaudin models was studied in [70J. 

Besides the approach proposed here, different recipes for quantizations of the Gaudin 
model were suggested, e.g., in [39] and [50]). The recipes proposed there are valid for 
the unreduced Gaudin models. Therefore, one needs to perform the quantum reduction 
procedure in order to relate the results of [19] and [50] with ours. 

• Higher spin chains. In this paper we consider the higher spin chains with orbits of 
the minimal rank at each site, i.e. S k £ slfc is conjugated to the element of the form 
diag(r, r, — (k — l)r). It is also interesting to describe the spectral duality for the 
generic GLj. spin chain. We hope to solve this problem in future publications. 

• Anisotropic chains. Another interesting generalization is induced by the five-dimen- 
sional AGT [7T] which implies a correspondence between the XXZ magnets (see |18j . 
[72j ) and a Gaudin- like model with relativistic (difference) dynamics. This latter would 
emerge, since on the conformal side one deals in this case with the q-Virasoro conformal 
block which implies a difference Schrodinger equation for the block with insertion of the 
degenerate field. 

A duality for the five dimensional quiver gauge theories was proposed in [18j . It relates 
the theories with the gauge groups SU{N) M ~ l and SU(M) N ~ 1 compactified on R 4 x S 1 , 
the radius of S 1 being i?5. At the level of integrable mechanics, the relevant system in 
this case is the XXZ spin chain. For the case of single SU(N) gauge group it can be 
written as follows [72] : 

w + Q 2 N(e 2<:/R5 )w- 1 = P N (e 2 ^ R5 ) , (7.9) 
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where Q 2 nW = Ilti( X ~ e 2m t/ R *){\ - e 2m i I s *) and P N {\) = Y\Zi( X ~ ea/Rs )) is 
a polynomial of degree N with the coefficients parameterizing the Coulomb branch of 
the vacuum moduli space. One observes that the curve has exactly the same form as 
in the four dimensional case, except that it is written in terms of the variables (w,e 2 ^). 
However, the SW differential is different, namely dS = (dw/w. 

Prom the results of [18] one gets the duality transformation for the spectral curve of the 
XXZ spin chain. It relates two different XXZ systems, corresponding to the two sides of 
the SU^Ny 1 - 1 o SU(M) N ~ 1 duality, that is, the iV-site g i M sp i n chain and the M-site 
gl N spin chain with the spin matrices of the minimal rank. More concretely, the duality 
exchanges the variables w e 2i */ Rb . The SW differential is manifestly invariant (up to a 
sign) under this transformation: dS = lne 2 ^ R5 dlnw = —d\ne 2< *^ R5 Inw. 

We are going to describe this duality explicitly in our future publication |74j . An extension 
to six dimensions (an elliptic extension of the differential operator in the Schrodinger 
equation versus the XYZ magnet) is also extremely interesting to construct. 
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